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Abstract: We investigate the possibility to measure the top quark mass using the collider 
variable rriT2 at the LHC experiment. Monte Carlo studies of are performed with the 
events corresponding to the dilepton decays of tt produced at the LHC with 10 fb^ 1 
integrated luminosity. Our analysis suggests that the top quark mass can be determined 
by the rriT2 variable alone with a good precision at the level of 1 GeV. 
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1. Introduction 

When the Large Hadron Collider (LHC) is turned on, it will serve as a 'top quark factory' 
(J, 0- The cross section for tt pair production at the LHC is estimated to be 833 pb at 
the NLO calculation ||, implying roughly 8 million it pairs per year at low luminosity run 
(10 /6 _1 /year). Such a large number of tt events will enable us to measure the top quark 
mass with high precision. 

Precision measurement of the top quark mass mt is desirable in many respects. For 
example, it would help to constrain the allowed Higgs boson mass in the Standard Model 
(SM). In general, it would affect the constraints on the allowed parameter space of var- 
ious models of new physics at the TeV scale, including the Minimal Supersymmetric 
Standard Model and technicolor-like models. The top quark mass measurement can be 
performed through various methods in different channels, which have their own advan- 
tage/disadvantage with different systematic uncertainties. In the overall, the accuracy of 
m t measured at the LHC is expected to be around 1 GeV Q . 

In the SM, top quark decays mostly into a b-quark and a W boson. The W boson 
then decays hadronically (W — > qq') or leptonically (W — > If). Depending on the W 
boson decay mode, the ti events are divided into three channels, i.e., the dilepton channel 
(both W bosons decay leptonically), the lepton plus jets channel (one W boson decays 
leptonically and the other hadronically) and the pure hadronic channel (both W bosons 
decay hadronically). 

The dilepton channel has a small branching fraction compared to the lepton plus jets 
channel and the pure hadronic channel. It also involves two missing neutrinos, which 
makes a direct event-by-event measurement of mt not possible. However, it has a cleaner 
environment, e.g. less combinatorial background and less jet energy scale dependence, 
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compared to other channels, therefore various approaches for an indirect measurement of 
mt with dilepton channel have been investigated Q. 

It has been shown that the collider variable m^2 H can be useful for the determination 
of new particle masses in the process in which new particles are pair produced at hadron 
collider and each of them decays into one invisible particle and one or more visible particles 
H> 0i §i §]• In this paper, we examine the possibility to determine the top quark mass 
using niT2 at the LHC experiment. For this, we perform three Monte Carlo studies of 
mj-2 for the process ti — > bl + vbl~v: the first which determines the endpoint value of the 
tot2 distribution for the neutrino mass m v = 0, the second to examine the functional 
dependence of m™| x on the trial neutrino mass m u 7^ 0, which would determine mt for a 
given value of the W boson mass mw-i an d the third which fits the rriT2 distribution to 
'template' distributions. Our analysis suggests that the top quark mass can be determined 
by the rriT2 variable alone with a good precision at the level of 1 GeV. 

In sec. 2, we briefly introduce the rriT2 variable for the dilepton decay of ti. The results 
of Monte Carlo studies are presented in sec. 3, and sec. 4 is the conclusion. 

2. Transverse mass and mri for top quark 

Let us consider a ti pair production and its subsequent decay at the LHC: 

pp -> tt -> bW + bW~ . (2.1) 

In case that one of the W bosons decays into leptons, one can consider the associated 
transverse mass of t — > blv, which is defined as 

m 2 T = m 2 +f7 ^ + 2 (£*f£^_ p «.p5,) j (2.2) 

where mu and denote the invariant mass and transverse momentum of the bl system, 
respectively, while m v and are the mass and transverse momentum of the missing 
neutrino, respectively. The transverse energies of the bl system and neutrino are defined 

as 

Et = \/\P b T\ 2 + m li and e t = y/\Pr\ 2 + m l- ( 2 - 3 ) 

If the other W boson decays into hadrons, i.e. for the process ti — > blubqq' , can 
be read off from the total missing transverse momentum p™ lss . One might then construct 
the rriT distribution of t — > blv from data, which can be used to determine the top quark 
mass mt as its shape and endpoint depend on mt- However, to determine p^. in the process 
ti — > blvbqq' ', one needs to measure the full final state momenta of i — ► bqq', which by itself 
would determine mt in event-by-event basis. At any rate, if one uses information from 
i — > bqq' to determine mt, the procedure involves more jets, which would result in larger 
uncertainties in the determined value of mt- 

A method to determine m t without using the hadronic decay of W is to construct mx2 
for the dilepton decay 

U = f(l)t(2) ^ 6 (l) Z (l) I/ (l) 6 (2) i (2) I/ (2) > (2 . 4) 
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Although each neutrino momentum cannot be measured in this case, still the total missing 
transverse momentum p™ lss = p^ 1 + p^ can be determined experimentally. The rriT2 
variable of each event is defined as 

mT2= vain max{m2 , } , (2-5) 

where m^) (i = 1,2) is the transverse mass of — > feWiWi/W, and the minimization is 
performed over the trial neutrino momenta p T constrained as 

Pt (1) +Pt 2) =Y>T iss - (2-6) 

The above definition of rriT2 indicates that rriT2 for m u = is bounded above by mt in 
the approximation ignoring the decay width of top quark. One might then determine rat 
as 

max/ n \ \ ( W bl W ( 2 ) bl ( 2 ) n\ 1 Co rrt 

m t = m T2 [m u = 0) = max m T 2(m^ , p T v % , p T , = 0) . (2.7) 

In fact, because of nonzero decay width, there can be certain amount of events which give 
rriT2 exceeding the physical top quark mass mt- Our Monte Carlo study suggests that such 
events do not spoil the sharp edge structure of the distribution with which one can 
determine mt rather precisely. Fig. [l] shows the top quark m,T2 distribution for m v = 



obtained from a parton level Monte Carlo simulation 1 using PYTHIA event generator [10] 
with an input top mass of mt = 170.9 GeV. One can see that mx2 tends to zero rapidly 
near the input top mass with a minor but long tail beyond the input mass which is mainly 
due to the nonzero top decay width 2 . 

One can consider the top quark mT2 defined as above for arbitrary trial neutrino mass 
which is not same as the true neutrino mass. In such case, mT2 is not only a function of 
the observable kinematic variables rm) and ' (i = 1,2), but also of the trial neutrino 
mass. Let m v denote the trial neutrino mass to distinguish it from the true neutrino mass 
m v = 0. The endpoint value of mx2 for generic m u , 

max/- s T / (1) bUl) (2) M(2) ~ \~\ fn o\ 

m T2 (m v ) = max mT2\m hl ' ,p T , m y bl ' ,p T w , m u ) , (2.8) 

appears to be a function of m„, and its functional form provides a relation between mt, 
the W boson mass my/, and the b quark mass m^. Using the result of Ref.||, one easily 
finds that m™2 X as a function of m v is given by 



where 

(»D 2 = m 2 b + l -{m 2 t -m 2 w - m\) + ^(m 2 - m 2 w - m 2 b Y - 4™^. (2.10) 



1 For simplicity, here we switched off the initial and final state radiations as well as the quark fragmen- 
tation process. 

2 Such a sharp edge structure of niT2 distribution at the input mass of the mother particle can be 
confirmed also in the tuti distribution for W + W~ — > l + vl~v. 
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This analytic expression of m™| ,x (?fijy) provides another way to determine mt, i.e. one can 
determine nit by fitting m™|' x (?fi l/ ) obtained from data to this analytic expression with the 
known values of mw and m&. 




Figure 1: 771^2 distribution obtained from partonic-level simulation. The input top quark mass of 
170.9 GeV is used for the simulation. One can find a sharp edge at the input top mass, with a small 
tail which is mainly due to the finite top quark decay width. 



3. Experimental feasibility 

Measuring the top mass using m.T2 in real experiment will suffer from a variety of uncer- 
tainty factors such as backgrounds, event selection cuts, finite jet energy resolution and 
combinatorial background. In order to check the feasibility of the rriT2 method at the LHC, 
we have generated Monte Carlo samples of tt events by PYTHIA [[[(]] with the CTEQ5L 
parton distribution function (PDF) The event sample corresponds to 10 fb -1 inte- 

grated luminosity. 

The generated events have been further processed with a modified version of fast 



detector simulation program PGS [12], which approximate an ATLAS or CMS-like detector 
with reasonable efficiencies and fake rates. The PGS program uses a cone algorithm for 
jet reconstruction, with default value of cone size AR = 0.5, where AR is a separation 
in the azimuthal angle and pseudorapidity plane. And the b-jet tagging efficiency e& is 
introduced as a function of the jet transverse energy and pseudorapidity, with a typical 
value of eb ~ 50% in the central region for high energy jets. 
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In the PGS, isolated leptons (electron and muon) are identified with some isolation 
cuts on the calorimeter activity around the lepton track [p^ j. For electrons, the isolation 
cuts are (i) ETISO / < 0.1, where ETISO is the total transverse calorimeter energy 
in a 3 x 3 grid around the electron candidate (excluding the candidate cell) and Et is 
the transverse energy of the electron candidate, (ii) PTISO < 5 GeV, where PTISO is 
the total pt of tracks (except the electron track) with pt > 0.5 GeV within a AR < 0.4 
cone around the electron candidate, and (iii) 0.5 < EP < 1.5, where EP is the ratio 
of the calorimeter cell energy to the px of the candidate track. For isolated muons, (i) 
PTISO < 5 GeV and (ii) ETRAT < 0.1125, where ETRAT is the ratio of E T in a 3 x 3 
calorimeter array around the muon (including the muon's cell) to the pt of the muon. 

The dilepton events are selected by requiring (A) only two isolated leptons of opposite 
charge with px > 25 GeV and \rj\ < 2.5, (B) dilepton invariant mass with \mu — mz\ > 5 
GeV, (C) large missing transverse energy E™ lss > 40 GeV, and (D) at least two b-tagged 
jets with pt > 30 GeV and \rj\ < 3.0. After this selection, 5133 events are survived among 
the 5.5 x 10 6 generated ti events (in which 1.8 x 10 5 are the dilepton events, considering 
only electrons and muons), leading to a selection efficiency of about 2.8% for the dilepton 
channel signal events. 

The main backgrounds might come from Z/~f* /W production with additional jets, 
diboson events with additional jets and bb events with misidentified leptons. We have 



generated the main background events using PYTHIA, ALPGEN [14( and AcerMC [15], 
and required the same selection cuts as the ti dilepton events. After the cuts, it turns 
out that those backgrounds are reduced to a negligible level. We will not include the 
background events in our further analysis, for simplicity. 

With two b-jets and two leptons in each selected event, there are two possible combina- 
tions for bl pairing. We calculated variable for each of the two possible bl combinations, 
and chose the smaller one as the final rriT2- This procedure closely follows the idea proposed 
inRef. g. 

Fig. |2] shows the resulting rriT2 distribution for the selected events. As anticipated, 
one can find an edge structure around mT2 = 170 GeV, on the distribution. We employ 
three methods to precisely determine the top quark mass from the rriT2 distribution, which 
will be discussed in the following three subsections. 

3.1 A fit near the end point 

Fig. ^ shows the rriT2 distribution obtained from the selected events for the neutrino 
mass m u = 0. It is fitted with an empirical function which consists of a linear function 
for signal distribution and an inverse linear function for background distribution. The fit 
range was chosen within ±O(10) bins around a plausible endpoint. Such fitting of the rriT2 
distribution results in 

m t = 171.1 ± 1.1 GeV, (3.1) 

which reproduces the input top quark mass of 170.9 GeV with a precision at the level of 1 
GeV. 
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Figure 2: rrvri distribution after event selection. The input value of top quark is m t = 170.9 GeV. 
A fit to the distribution near endpoint region is also shown, providing a fit value of m t = 171.1 ± 1.1 
GeV. 

To estimate possible systematic error associated with the fitting procedure, we have 
repeated the fitting with two linear functions for both signal and background distributions. 
The resulting top quark mass is then given by mt = 169.9 ± 1.8 GeV, showing a mass shift 
of 1.2 GeV. Systematic error from the fitting procedure might be improved by considering 



a template binned likelihood fit, which will be discussed in subsection |3,3 . 

Absolute jet energy scale also affects the determination of the top mass. The b-jet 
energy scale is assumed to be known within 1% accuracy. It is found that the 1% variation 
of the jet scale leads to a shift of the resulting top mass of 0.5 GeV. 

Uncertainty due to initial state radiation (ISR) is estimated by comparing the nominal 
data (with ISR switched on) to the one which is generated while switching off ISR. The 20% 
of the resulting top mass shift is found to be 0.4 GeV, which is taken as the systematic error 
from ISR uncertainty Q . The same approach to final state radiation induces a systematic 
error of 0.7 GeV. 

For systematic error from PDF uncertainty, it is found that the use of CTEQ3L 
(GRV94L) PDF, instead of the default CTEQ5L PDF, leads to a shift of the central top 
mass of 0.3 (1.3) GeV, with a suitable choice of fit range. 

3.2 Endpoint as a function of trial neutrino mass 

As we have discussed in section ^, the endpoint of rriT2 distribution can be considered as 
a function of a trial neutrino mass, if we use a trial neutrino mass m v ^ for the mr2 



- 6 - 




Figure 3: (a) An example of '111x2 distribution with a trial neutrino mass. Here, the trial mass 
is set to rh v = 80 GeV. (b) The maximum of rriT2 as a function of trial neutrino mass fh v . Also 
shown is the fit of the data points to theoretical curve ( 2. 7) considering m t as a free parameter. 

calculation. Using the selected dilepton decays of ti, we constructed the rriT2 distributions 
for different choices of fh u . Fig. |3|(a) shows the mT2 distribution for m v = 80 GeV. Here 
we also performed a fit to the m^2 distribution with a linear function for signal and an 
inverse linear function for background. The maximum of ntTi is then determined to be 
m T2 X = 232.6 ± 1.5 GeV for m v = 80 GeV. The m™| x as a function of m u is shown in 



Fig. |3|(b). Fitting the data points to the theoretical curve (2J3) considering mt as a free 



parameter while using mw = 80.45 GeV and to?, = 4.7 GeV, we obtain 

mt = 170.5 ± 0.5 GeV, (3.2) 

which is quite close to the input top quark mass mt = 170.9 GeV. The uncertainty due to 
a variation of m^ is negligible as it is of Oim^bm^jmt). To check the effect of the W boson 
mass, we repeated the fitting procedure while varying mw by ±0.5 GeV. The resulting 
shift of mt turns out to be negligible. 

3.3 Template binned likelihood fit 

Perhaps the most reliable way to determine mt using mT2 is to employ the template 
binned likelihood fit. For this, we attempted to fit the mT2 distribution of the 'nominal 
data' (which was generated with mt = 170.9 GeV) to 'templates'. Here, a template means 
a simulated mx2 distribution with an input top quark mass different from 170.9 GeV. The 
templates were generated with input top quark mass between 166 GeV and 176 GeV, in 
steps of 1 or 0.5 GeV, using the same PYTHIA+PGS Monte Carlo programs as the case 
of nominal data sample. 
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Figure 4: (a) Three representative tuti distributions for the nominal data (points) and two tem- 
plates with mt = 166 GeV (blue solid) and mt = 176 GeV (red solid), respectively, (b) The negative 
logarithm of the likelihood ratio C/C^^ as a function of m t for the TOT2 fit. 

Fig. H|(a) shows three representative rriT2 distributions for the nominal data (points) 
and two templates with mt = 166 GeV (blue solid) and mt = 176 GeV (red solid), re- 
spectively. Each template distribution is normalized to make the total number of events 
is same as that of the nominal data. One can notice that those three mx2 distributions 
are well separated from each other, showing the sensitivity of the distribution to the 
input top quark mass. 

Each template distribution is compared to the nominal data distribution for a calcu- 
lation of the logarithm of the binned likelihood. The binned likelihood is defined as the 
product of the Poisson probability for each bin over the N bins in the fit range: 



N 



-rriiTii 



i=l 



where rij and m; are the event numbers at the i-th bin in the distributions of the nominal 
data and the normalized template, respectively. The minimum of —InC gives the best fit 
value of the top quark mass. We have chosen the la deviated value of the top quark mass 
as the one increasing — ln£ by 1/2. 

We fit the mT2 distribution of nominal data to templates in the range 100 GeV < 
mx2 < 180 GeV. The result of the likelihood fit for m?2 distributions is shown in Fig. 
H](b), where the negative logarithm of the likelihood ratio £/£ max as a function of mt is 
depicted. The £ max is the maximum likelihood which was determined as the minimum of 
a parabola fit to the — ln£ distribution. The top quark mass resulting from our template 
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likelihood fit is given by 



m t = 170.3 ± 0.3 GeV, (3.4) 

which reproduces well the input top quark mass with a small statistical error. 

Although a detailed analysis of systematic uncertainties in the template fit method is 
beyond the scope of this work, we expect that systematic errors from b-jet energy scale, 
ISR/FSR and PDF are also at the level of 1 GeV as those in the endpoint fit method 



discussed in subsection 3.1 



4. Conclusion 



We have examined the possibility to determine the top quark mass using the m,T2 distribu- 
tion of the dileptonic decay channel of ti events at the LHC. For this, we have performed 
three Monte Carlo studies for the events produced at the LHC with 10 fb^ 1 integrated 
luminosity: the first to fit the rriT2 distribution near the end point (for the neutrino mass 
m v = 0) with an empirical function, the second to fit the functional dependence of m™| x 
on the trial neutrino mass rh u ^ 0, and the third to perform a template binned likelihood 
fitting. It is found that the top quark mass can be determined by the mr2 variable alone 
with a good precision at the level of 1 GeV. 
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